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Abstract 

Recently it was pointed out that for the evaluation of the numerically dominant 
pion-exchange contribution to the hadronic light-by-light scattering correction in 
the muon g — 2, a fully off-shell pion-photon-photon form factor should be used. 
Following this proposal, we first derive a new short-distance constraint on the off- 
shell form factor which enters at the external vertex for the muon g — 2 and show 
that it is related to the quark condensate magnetic susceptibility in QCD. We then 
evaluate the pion-exchange contribution in the framework of large- Nc QCD using an 
off-shell form factor which fulfills all short-distance constraints. With a value for the 
magnetic susceptibility as estimated in the same large- Ac framework, we obtain the 
result <v byL; *" = (72±12)xl0- n . Updating our earlier results for the contributions 
from the exchanges of the n and rf using simple vector-meson dominance form 
factors, we obtain a]l hyL ' PS = (99±16) x 10" 11 for the sum of all light pseudoscalars. 
Combined with available evaluations for the other contributions to hadronic light- 
by-light scattering this leads to the new estimate a L b y L > had _ ^16 ± 40) x 10~ n . 
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1 Introduction 



The hadronic light-by-light scattering contribution to the muon g — 2 has a long and 
troubled history. The relevant physics involves the nonperturbative regime of QCD below 
about 1 — 2 GeV. Furthermore, no direct experimental information is available, in contrast 
to the hadronic vacuum polarization contribution to the g — 2, which can be related to 
the cross section e + e~ — > hadrons. Therefore various models have been used over the 
years, starting first with a simple loop of some constituent quarks [1] . Later more realistic 
hadronic models with exchanges of light pseudoscalars, scalars and other resonances and 
loops with charged pions have been employed [2]. However, the coupling of hadrons 
to photons will, in general, involve some form factors which are very model-dependent 
[p — 7 mixing as in vector- meson dominance ( VMD) models] . In the absence of any direct 
experimental checks, the size and even the sign of the light-by-light scattering contribution 
to the muon g — 2 was therefore uncertain for a long time. Actually, the sign has changed 
several times over the years due to some errors in the complicated calculations. 

In Ref. [3] a systematic approach was proposed, based on the chiral expansion [I] and 
the large- Nq counting [S] of the various contributing diagrams. Soon afterwards, two 
very extensive evaluations appeared, Refs. P, [TJ [8] , based on slightly different hadronic 
models. However, they both had a sign error in the numerically dominating pseudoscalar- 
exchange contribution as was pointed out a few years later in Refs. [H HD] and confirmed 
in Refs. [IB CE2J H3] ■ 

Reference [5] mainly concentrated on the neutral pion-exchange contribution where the 
pion-photon-photon form factor jF^o*^. enters. In general, low-energy hadronic models 
for form factors, e.g. based on some constituent quark model or on some resonance La- 
grangian, do not satisfy all the large momentum asymptotics required by QCD. Using 
these form factors in loop diagrams thus leads to cutoff-dependent results. Even if the 
cutoff is varied in a reasonable range, e.g. ~ 1 — 2 GeV, the corresponding model un- 
certainty is completely uncontrollable. In order to eliminate (or at least reduce) this 
cutoff dependence, new models for J 7 ^ *^* were proposed in Ref. [Hj and then applied 
to hadronic light-by-light scattering in Ref. |9j. These models are based on the large- Nc 
picture of QCD [5], where, in leading order in Nc, an (infinite) tower of narrow resonances 
contributes in each channel of a particular Green's function. The low-energy and short- 
distance behavior of these Green's functions is then matched with results from QCD, 
using chiral perturbation theory [I] and the operator product expansion (OPE) [15], re- 
spectively. Based on the experience gained in many examples of low-energy hadronic 
physics, and from the use of dispersion relations and spectral representations for two- 
point functions, it is then assumed that with a minimal number of resonances in a given 
channel one can get a reasonable good description of the QCD Green's function in the 
real world (minimal hadronic Ansatz) . Often only the lowest-lying resonance is considered 
[lowest-meson dominance (LMD)] [THl CEO US UHl [20], as a generalization of VMD. 

Reference [9] obtained the result a^ byL;7r ° = (58±10) x KT 11 for the pion and aJ; byL;PS = 
(83± 12) x 10 -11 for the sum of all light pseudoscalars tt°, rj, and rj'. These results are close 
to the (sign corrected) values a^ byL ^° = (59 ± 9) x KT 11 [aj; byL;PS = (85 ± 13) x KT 11 ] 



1 



obtained in Ref. [B] and a^ byL; ^ = (57 ± 4) x 1CT 11 [a pbyL;PS = (82.7 ± 6.4) x 1CT 11 ] in 
Refs. [HIE]- The results for the (corrected) full contributions at that time read a L b y L ; had = 
(83 ± 32) x 10" 11 [6] and a ;; byL;had = (89.6 ± 15.4) x 10" 11 [HE]. 

Later Ref. [21] pointed out that some additional QCD short-distance constraint was 
not taken into account for the exchanges of pseudoscalars and axial- vector resonances in 
Refs. [HI 13 [HI [5j. The authors of Ref. [21] argued that if one imposes this constraint, no 
momentum-dependent form factor can be present at the external vertex which couples 
to the soft photon relevant for the magnetic moment. In the absence of such a form 
factor, Ref. [21] then got enhanced results compared to the earlier evaluations a^ byL;7r = 
(77 ± 7) x KT 11 [aJ; byL;PS = (114 ± 10) x 1CT 11 ] and a |; byL;had = (136 ± 25) x 10 _u . As 
discussed in Ref. [22] , a part of the enhancement of the result for a L b y L ; had - n j^ e f j^jj j s 
actually due to a different treatment of the axial vectors (ideal mixing instead of nonet 
symmetry) and the omission of the negative contributions from scalar exchanges and the 
charged pion loop. It is therefore not entirely related to the new short-distance constraint. 
Thus the slightly lower estimate a^ byL;had = (110 ± 40) x 10~ n has been employed in the 
reviews [221 [23]. Very recently, the value a L b y L ; had _ (i05±26) x 10~ n has been proposed 
in Ref. [21]. 

However, recently Refs. [251 [2H] stressed the fact that one should actually use fully 
off-shell form factors for the evaluation of the light-by-light scattering contribution. This 
seems to have been overlooked in the recent literature, in particular, in Refs. P, [271 EH 
[22]. The on-shell form factors as used in Refs. [91 [27] actually violate four-momentum 
conservation at the external vertex. While Ref. [21] had already pointed out this violation 
of momentum conservation at the external vertex, they then only considered on-shell pion 
form factors, an approximation which yields the so-called pion-pole contribution and not 
the more general pion- exchange contribution with off-shell form factors. Putting the pion 
on-shell at the external vertex automatically leads to a constant form factor. 

In the present paper we revisit the pion-exchange contribution in view of the obser- 
vations made in Refs. [251 [2S]- We first derive a new QCD short-distance constraint for 
the off-shell form factor which enters at the external vertex and show that it is related 
to the quark condensate magnetic susceptibility in QCD. We also comment on how our 
short- distance constraint is connected with the one derived in Ref. |21j. In the second 
part we evaluate the pion-exchange contribution in the framework of large- Nc QCD with 
off-shell form factors both at the internal and the external vertex, taking into account the 
new short-distance constraint and an estimate for the magnetic susceptibility in QCD in 
the same large- Nc framework. 

Strictly speaking, the identification of the pion-exchange contribution is only possible, 
if the pion is on-shell (or nearly on-shell). If one is (far) off the mass shell of the exchanged 
particle, it is not possible to separate different contributions to the g — 2, unless one 
uses some particular model where for instance elementary pions can propagate. In this 
sense, only the pion-pole contribution with on-shell form factors can be defined, at least 
in principle, in a model-independent way, although the numerical result will in general 
still depend on the model used for the on-shell form factors, unless one would know the 
"true" form factors. On the other hand, the pion-pole contribution is only a part of the 
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full result, since in general, e.g. using some resonance Lagrangian, the form factors will 
enter the calculation with off-shell momenta. In this respect, we view our evaluation as 
being a part of a full calculation of hadronic light-by-light scattering using a resonance 
Lagrangian whose coefficients are tuned in such a way as to systematically reproduce the 
relevant QCD short- distance constraints, e.g. along the lines of the resonance chiral theory 
developed in Ref. [2"g] . 

We should mention that recently another paper appeared [2H] which evaluates the 
pion-exchange contribution using an off-shell form factor based on the nonlocal chiral 
quark model, obtaining the result aJ^ byL;7r = (65 ± 2) x 10~ n . We will comment on that 
paper below. 

This paper is organized as follows. Section [2] contains the starting point for the cal- 
culation of the pion-exchange contribution to the muon g — 2, including the definition of 
the pion-photon-photon form factor .F 7r o* 7 » T *- We also discuss the issue of using on-shell 
or off-shell form factors. In Sec. [3] we discuss several experimental and theoretical con- 
straints on the form factor. In particular, we derive a new short-distance constraint on 
the off-shell form factor at the external vertex in hadronic light-by-light scattering. In 
Sec. H] we present a new evaluation of the pion-exchange contribution in the framework of 
large- Nq QCD and give some updated estimates for the r\ and rj exchange contributions 
using simple VMD form factors. We end with discussions and conclusions in Sec. El In 
the appendix we give a parametrization of the numerical result for the pion-exchange 
contribution for arbitrary parameters of our model for the off-shell form factor. 

2 The pseudoscalar-exchange contribution 

The numerically dominating contributions to hadronic light-by-light scattering are due to 
the neutral pseudoscalar-exchange diagrams shown in Fig. [TJ 




Figure 1: The pseudoscalar-exchange contributions to hadronic light-by-light scattering. 
The shaded blobs represent the off-shell form factor JF PS * 7 . 7 * where PS = 7r°, rj, r}'. 

We first concentrate on the exchange of the neutral pion. The key object which 
enters the Feynman diagrams in Fig. [1] is the off-shell pion-photon-photon form factor 
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■7^.o* 7 * 7 » + q 2 ) 2 , qf, q 2 ) which is defined, up to small mixing effects with the states r\ 
and rj', via the Green's function (VVP) in QCD 



d 4 xd 4 ye^ x+ ^yH0\T{j^x)jMP 3 m\0) 



^ (gi + 92) 2 - ml 



(i) 



Here = [ijjQ^^ijj)^) [ip = (u, d,s), Q = diag(2, —1, — 1)/3 the charge matrix] is the 

light quark part of the electromagnetic current and P 3 = ipi^Y^ = (uijsu — dijsd) /2. 
Note that we denote by (ipijj) the single flavor bilinear quark condensate. The form factor 
is of course Bose symmetric ^ r 7r o* 7 * 7 »((g 1 + q 2 ) 2 , q\, q 2 ) = J r 7r o* 7 * 7 *((gi + g 2 ) 2 , ?f> as the 
two photons are indistinguishable. 

The corresponding contribution to the muon g — 2 may be worked out with the re- 
sult [90 

<i 4 gi d 4 g 2 1 

(2vr)4 (2tt)4 g 2 g 2 2 ( gi + q 2 f[{p + gi ) 2 - ™ 2 ][(P - <? 2 ) 2 - m 2 ] 

2i(gi,g 2 ;p) 



LbyL;7r° 



X 



+ 



^V>Vt* (?2> 9i> (ft + ^) 2 ) ^V0V7(92 J ?2> °) 



?| -m 2 



Jv>Vt*((?i + 9 2 ) 2 ' 9i> 92) ^> 7 *t((9i + 92) 2 , (qi + g 2 ) 2 , 0) 



2iyi 



92; p) 



(2) 



with 



Ti(qi,q 2 ;p) 



16 



(p- 9i) (P* 92) (9i • 92) 



16 



(P • 92) 2 ql 



16 



g (P ■ Qi) (Qi ■ 92) ql + 8(p • q 2 ) q\q\ - y (p • 92) (?i • 9 2 ) 2 



T%{q\,qi\V) 



1 xu 2 2 2 iu 2 / \2 

+ y m M 9i9 2 - -fe) 

16 16 



(3) 



y (p ' 9i ) (P • 92) (<?i -92) - y (P • <?i) 2 9 2 
+ 3 (P • 9i) (9i • 92) <?2 + g (P • 9i) 9i 9 2 

8992 2 / \2 

+ 3^9l?2 - g"M9l -92) , 



(4) 



where p 2 = m 2 and the external photon has now zero four-momentum. The first and 
the second graphs in Fig. [1] give rise to identical contributions, leading to the term with 
Ti, whereas the third graph gives the contribution involving T 2 . The factor T 2 has been 
symmetrized with respect to the exchange q\ <-> —q 2 . 



1 To be precise, the corresponding expression with on-shell form factors is given in Ref. [S]. 
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Instead of the expressions in Eq. Refs. [HI EZ] and maybe also earlier works, 
considered on-shell form factors, e.g. for the term involving T 2 one would write [26J 

^VtK>?i.?2 2 ) x -TvVtC 771 *' (ft + ^) 2 ,0). (5) 

Often the first argument of the on-shell form factor is omitted in the literature, i.e. the 
form factor is written as ^ r 7r o 7 * 7 * (g 2 , g|) = ^ r 7r o 7 * 7 »(m^, gf , gf). Although pole dominance 
might be expected to give a reasonable approximation, it is not correct as it was used 
in those references, as stressed in Refs. [2"Tj |2"B"1 The point is that the form factor 
sitting at the external photon vertex in the pole approximation jF 7r o 7 * 7 (m^., (gi + g 2 ) 2 ,0) 
for (gi + g 2 ) 2 7^ ^ 2 violates four- momentum conservation, since the momentum of the 
external (soft) photon vanishes. The latter requires JF^o* 7 * 7 ((gi + g 2 ) 2 , (li + <?2) 2 ,0). In 
order to avoid this inconsistency, Ref. [2T] proposed to use instead 

^V7*K>9i'% 2 ) x -^ 7 (™ 2 ,™ 2 ,0) , (6) 

i.e. a constant form factor at the external vertex, which is given by the Wess-Zumino- 
Witten (WZW) anomaly [30]. The absence of a form factor at the external vertex in the 
pion-pole approximation follows automatically, if one carefully considers the momentum 
dependence of the form factor. This procedure is also consistent with any quantum field 
theoretical framework for hadronic light-by-light scattering, for instance, if one uses a 
(resonance) Lagrangian to derive the form factors, and where a different treatment of the 
internal and external vertex, apart from the kinematics, is not possible. On the other 
hand, taking the diagram more literally, would require 

F 7r o* rr ((q l + q 2 ) 2 ,ql,ql) x ^^((gi + g 2 ) 2 , (<?i + ?2) 2 , 0) , (7) 

as the more appropriate amplitude, see Eq. (j2J). References [251 EE] advocate the use of 
fully off-shell form factors at both vertices and we will follow this approach in the rest of 
this paper. The difference to the procedure adopted in Ref. [21] will be important when 
we discuss their short- distance constraint. 



3 Experimental and theoretical constraints on the 
pion-photon-photon form factor 

The form factor ^vo* 7 * 7 *((gi + g 2 ) 2 , g 2 , gf) defined in Eq. (TJ]) is determined by nonpertur- 
bative physics of QCD and cannot (yet) be calculated from first principles. Therefore, 
various hadronic models have been used in the literature, sometimes combined with short- 
distance constraints from perturbative QCD at high momenta. At low energies, the form 
factor is normalized by the decay amplitude, A(n° — > 77) = e 2 jF 7r o 77 (m 2 r , 0, 0) in our 
conventions. In the chiral limit, m q — > 0, q = u, d, s, this amplitude is fixed by the WZW 
anomaly 

^»> ( ^ 77) = __0|. (8) 
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For massive light quarks, this expression receives corrections. In particular, the pion decay 
constant in the chiral limit, F , is replaced by its physical counterpart F n = F [l + O (m q )\. 

^°-77) = - I ^|r[l + 0K)]. (9) 

It turns out that the additional O (m q ) corrections in this relation are numerically 
small [31], so that one may drop them to a good approximation. The measured de- 
cay width r(7T° -> 77) = (7.74 ± 0.6) eV [32] is then well reproduced for F w = 92.4 MeV. 
Therefore, all hadronic models for the form factor have to satisfy the low-energy constraint 

jr O77 K )0 ,0) = -^^. (10) 

Sometimes this normalizing value is used to define a constant "WZW form factor" 
^JaSJL, ( (qi + (fc) 2 , qf, q£) = —Nc/{127r 2 F n ). This notion of a constant form factor is, 
however, very misleading. For off-shell momenta away from the physical point in Eq. (fTOl 
the value of this WZW form factor has no physical meaning. Recall that the WZW effec- 
tive Lagrangian only yields the first term in the low-energy and chiral expansion of the 
corresponding (VVP) Green's function. 

For an on-shell pion, there is also experimental data available for one on-shell and one 
off-shell photon, from the process e + e~ — *> e + e~7r°. Several experiments [231 El] thereby 
fairly well confirm the Brodsky-Lepage [35] behavior for large Euclidean momentum 

2F 

lim j;o 7 , 7 (m 2 ,-Q 2 ,0) ~ - — f (11) 

and any satisfactory model should reproduce this behavior. 

Apart from these experimental constraints, any consistent hadronic model for the 
off-shell form factor ^ r 7r o* 7 » 7 *((gi + q2) 2 , qf, q 2 ) should match at large momentum with 
short-distance constraints from QCD that can be calculated using the OPE. In Ref. [H] 
the short- distance properties for the three-point function (VVP) in Eq. ([T]) in the chiral 
limit and assuming octet symmetry have been worked out in detail (see also Refs. [T7] I2TJ] 
for earlier partial results). At least for the pion the chiral limit should be a not too bad 
approximation^; however, for the rj and, in particular, for the non-Goldstone boson rj 
further analysis will be necessary. 

It is important to notice that the Green's function (VVP) is an order parameter of 
chiral symmetry. Therefore, it vanishes to all orders in perturbative QCD in the chiral 
limit, so that the behavior at short distances is smoother than expected from naive power 
counting arguments. Two limits are of interest. In the first case, the two momenta become 
simultaneously large, which in position space describes the situation where the space-time 
arguments of all three operators tend towards the same point at the same rate. To leading 

2 As pointed out in Ref. [36], the integrals in Eq. ([2|) are infrared safe for — > 0. This can also be 
seen within the effective field theory approach to light- by- light scattering proposed in Refs. pTCTl [T5] . 
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order and up to corrections of order O (a s ) one obtains the following behavior for the form 
factoiB 



lim ^o Vr ((Aft + M 2 )\ (A 9l ) 2 , (Ag 2 ) 2 ) = f ^ + gi + 2 ( f + 92)2 + O (± ) . (121 

A^oo OA Q1Q2 \^ 



The second situation of interest corresponds to the case where the relative distance 
between only two of the three operators in (VVP) becomes small. It so happens that 
the corresponding behaviors in momentum space involve, apart from the correlator (AP) 
which, in the chiral limit, is saturated by the single-pion intermediate state 

d i xe vpx (Q\T{A a (x)P b ^)}\Q) = 5 ab ^) ^ (13) 

(we denote by (ipip)o the single flavor bilinear quark condensate in the chiral limit), the 
two-point function (VT) of the vector current and the antisymmetric tensor density 



5 ab (n V TW(p) = / d 4 xe^(0\T{VVx)^a pcT -^(0)}\0) , (14) 



with a ^ = §[7(0,70-] (the similar correlator between the axial current and the tensor 
density vanishes as a consequence of invariance under charge conjugation). Conservation 
of the vector current and invariance under parity then give 

(n V T) MP a(p) = (p p v^ - PaVw) n V T(p 2 ) • (15) 

When the space-time arguments of the two vector currents in (VVP) approach each 
other, the leading term in the OPE leads to the Green's function (AP) and the short- 
distance behavior of the form factor reads 

IF 11 / 1 \ 

Km Jv» Vr (ql (A gi ) 2 , (q 2 - \ gi ) 2 ) = + O i-^j . (16) 

Further important information on the on-shell pion form factor has been obtained in 
Ref. [37J based on higher-twist terms in the OPE and worked out in [38]. In the chiral 
limit one obtains the behavior 

J- 7r o rr (0,(Ag 1 ) 2 ,(Ag 1 ) 2 ) 8_ 2r?2 f 1 ,8 5 2 ^(1 



& 7 >^(o,o,o) = ~r' F « { av + » av + Uj ) • < 17 > 

where S 2 parametrizes the relevant higher- twist matrix element. The sum rule estimate 
performed in [SB] yields the value 5 2 = (0.2 ± 0.02) GeV 2 . 



3 In the chiral limit, the relation between the off-shell form factor and the single invariant func- 
tion Hv which appears in (VVP) is given by ^ r o« 7 » 7 ,((gi + q 2 ) 2 , ql, q%) = — (2/3)(i i o/ '(^)o)(qi + 
q2) 2 'Hv(ql, ql, (qi + ffii) 2 ); see Ref. [H] for details. 
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On the other hand, when the space-time argument of one of the vector currents in 
(VVP) approaches the argument of the pseudoscalar density one obtains the relation 



lim Jv, Vr ((Agi + q 2 ) 2 , (Xqi) 2 ,ql) = -^^^n VT (g 2 2 ) + O ( - ) . (IS) 



A— »oo 



3(^)o yz ' VA 



In particular, at the external vertex in light-by- light scattering in Eq. (J2J), the following 
limit is relevant 



— l(^ nvT(o)+o G 



lim J^WtAft) 2 , (A 9l ) ! ,0) = --^^nvrtO) +0 ( T ) . (19) 



Note that there is no falloff in this limit, unless Hvt(O) vanishes. 

As pointed out in Ref. [39], the value of IE/tG^ 2 ) at zero momentum is related to the 
quark condensate magnetic susceptibility x m QCD in the presence of a constant external 
electromagnetic field, introduced in Ref. 



(0\qa^q\0) F = ee qX (^)o F^, (20) 

with e u = 2/3 and = —1/3. With our definition of Hvt in Eq. (fT4l) one then obtains 
the relation (see also Ref. [4T] ) 

n VT (o) = — x, (21) 

and therefore the behavior at the external vertex from Eq. (fl~9l) can be rewritten as 

Km (Agx) 2 ,0) = ^ x + <5 (22) 

Unfortunately there is no agreement in the literature what the actual value of \ should 
be. In comparing different results one has to keep in mind that \ actually depends on the 
renormalization scale fx. In Ref. [ID] the estimate x(a* = 0-5 GeV) = — (8. 16jlf ) GeV -2 
was given in a QCD sum rule evaluation of nucleon magnetic moments. This value was 
confirmed by the recent reanalysis [12] which yields x = — (8.5 ± 1.0) GeV -2 , although 
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no scale fx has been specified. A similar value x — ~~ ^ r c/(4vr 2 F^) = —8.9 GeV 
obtained in Ref. [IB]- From the explicit expression of x it is n °t immediately clear what 
should be the relevant scale fx. Since pion dominance was used in the matching with the 
OPE below some higher states, it was argued in Ref. |43j that the normalization point is 
probably rather low, fx ~ 0.5 GeV. Calculations within the instanton liquid model yield 
X ILM (/i ~ 0.5 — 0.6 GeV) = —4.32 GeV -2 [H], where the scale is set by the inverse average 
instanton size p~ l . The value of xi^^o — 42 MeV at the same scale obtained in Ref. [B] 
agrees roughly with the result 35 — 40 MeV from Ref. [15] derived in the same model. 

The leading short-distance behavior of the two-point function IIvt i n Eq. (fT5"]) is given 
by [H] (see also Ref. [IS]) 



1 ^ttk + o (1 

A^o Vivv '' r/ > \2 p 2 \\ 



lim E VT ((A P ) 2 ) = -±^m + ol^) . (23) 
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Assuming that IIvt(p 2 ) is well described by the multiplet of the lowest-lying vector mesons 
(LMD) and satisfies the OPE constraint from Eq. QEfy leads to the Ansatz [4TI 1591 [T4] 

n v W) = -(^)o-A^- ( 24 ) 

v lvl v 

Using Eq. (ED then leads to the estimate x LMD = -2/M$ = ~ 3 - 3 GeV -2 [47]. Again, it 
is not obvious at which scale this relation holds. In analogy to estimates of low-energy 
constants in chiral Lagrangians [28], it might be at p = My, although in principle the 
renormalization scale of x is n °t related to the one of low-energy constants; see the 
discussion in Ref. [IB]. This LMD estimate was soon afterwards improved by taking into 
account higher resonance states (p', p") in the framework of QCD sum rules, with the 
results x(0.5 GeV) = -(5.7 ± 0.6) GeV~ 2 [39J and x(l GeV) = -(4.4 ± 0.4) GeV" 2 [49]. 
A more recent analysis [50] yields, however, a smaller absolute value %(1 GeV) = — (3. 15± 
0.30) GeV" 2 , close to the original LMD estimate. Further arguments for the latter value 
are also given in Ref. [4T] and references therein, by studying the coupling of the tensor 
current to the p meson. For a quantitative comparison of all these estimates for x we 
would have to run them to a common scale, for instance, 1 GeV, which can obviously 
not be done within perturbation theory starting from such low scales as p = 0.5 GeVjj 
Finally, even if the renormalization-group running could be performed nonperturbatively, 
it is not clear what would be the relevant scale p in the context of hadronic light-by-light 
scattering. 

A short-distance constraint on the pion-exchange contribution to the hadronic light- 
by-light scattering correction in the muon g — 2 itself was derived in Ref. [21]. The relevant 
kinematical configuration for the s-channel exchange of the pion is shown in Fig. [21 In 




Figure 2: Pion exchange in the s channel in hadronic light-by-light scattering. The photon 
with zero momentum g 4 = represents the external soft photon for the corresponding 
contribution to the muon g — 2. 

general one has q\ + q 2 + q% + q± = 0, but for the muon g — 2 the soft photon limit 

A further complication arises in comparisons with papers from the early 1980's because not only 
(i = 0.5 GeV was frequently used, but also 1-loop running with a low Aq^, 3 = 100 — 150 MeV, whereas 

more recent estimates yield A^-T 3 = 346 MeV (at 4- loop) [5"T] . 
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q4 — > will be relevant. The authors of Ref. then consider the limit qf ~ q\ ^ lii 
where q% = —{q\ + (72). Since in this limit the leading term in the OPE of the two vector 
currents associated with the momenta q\ and q2 yields the axial- vector current, they 
can relate the matrix element (VVVI7) which enters for the muon g — 2 to the famous 
anomalous triangle diagram (AVI7) [52], which is highly constrained; see Refs. [531 H3]. 
From this they deduce that no momentum-dependent form factor should be used at the 
external vertex, but only a constant factor. They thus obtain the following intermediate 
expression for the light-by-light scattering amplitude^: 

A n o = JL ^Vrfegl) if ^J n{ f pa f°P) + permutations, (25) 
Zr w q 3 m n 

where ft 1 " = q%e\ — g^ef and fi ;fJiL/ = \e^ vpa f pa for i = 1,2,3 denote the field strength 
tensors of the internal photons with polarization vectors The field strength tensor of 
the external soft photon is defined similarly by f pu = q p e u A — q\e\. Except in f pa the limit 
g 4 — > is understood in Eq. ([2"5lh in particular, in f^ p and in the pion propagator. 

Note the absence of a second form factor ^ r 7r o 7 * 7 (g|, 0) in Eq. f[2"5l) at the external vertex. 
The authors of Ref. [21] rightly point out that such a momentum-dependent form factor 
at the external vertex would violate momentum conservation and criticize the procedure 
adopted in earlier works [61 [H El [9] . However, it is obvious from their expressions [Eq. (18) 
in Ref. [2Tj]. reproduced here in Eq. (I25p . that they only consider the on-shell pion form 
factor F^o^j*^, q%) = J r 7r o 7 * 7 « (m^., qf, q%) at the internal vertex and not the off-shell pion 
form factor •F 7r o* 7 . 7 *(q , 3, qf, ql). Note that the expression in Eq. (l25l) has to be compared 
with the term involving T2 in Eq. (T5]). Therefore, contrary to the claim in Ref. [21], they 
only consider the pion-pole contribution to hadronic light-by-light scattering and not the 
pion-exchange contribution which involves fully off-shell form factors at the internal and 
the external vertex. Actually, also a second argument in Ref. [21J [after Eq. (20) there] in 
favor of a constant form factor at the external vertex is clearly based on the use of on-shell 
form factors. The use of a nonconstant on-shell form factor ^ r 7r o 7 * 7 (q , |, 0) at the external 
vertex would lead, together with the pion propagator, to an overall l/g| behavior, since 
jF 7r o 7 * 7 (g|, 0) ~ I/?!, f° r l ar g e ih according to Brodsky-Lepage; see Eq. (fTTT) . This would 
contradict the l/g| behavior observed in Eq. (1251) (apart from f^ p ). 

4 New evaluation of the pseudoscalar-exchange con- 
tribution in large- Nq QCD 

In the spirit of the minimal hadronic Ansatz for Green's functions in large-A^c QCD, 
on-shell jF 7r o 7 * 7 * (m^, qf, q\ ) and off-shell form factors ^ r 7r o* 7 * 7 *((gi + q 2 ) 2 , qf, q%) have been 
constructed in Ref. [14J. They contain either the lowest-lying multiplet of vector reso- 
nances (LMD) or two multiplets, the p and the p' (LMD+V). Both Ansatze fulfill all 



5 We have rescaled the form factor in Eq. (18) in Ref. [3T] to agree with our normalization in Eq. (jTUJ) 
and used Minkowski space notation. 
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the OPE constraints from Eqs. (1121) . (I16p and (1181) ; however, the LMD Ansatz does 
not reproduce the Brodsky-Lepage behavior from Eq. (fiT|) . Instead it behaves like 
J 7 ^, ' {ml, -Q 2 ,0) ~ const. The l/Q 2 falloff can be achieved with the LMD+V Ansatz 
with a certain choice of the free parameters; see below. Note that it might not always 
be possible to satisfy all short- distance constraints, in particular from the high-energy 
behavior of form factors, if only a finite number of resonances is included; see Ref. [53] . 
The on-shell form factors were later used in Ref. [9] to evaluate the pion-pole contribution; 
see also Ref. [2"T] . 

In the following, we reevaluate the pion-exchange contribution using off-shell LMD+V 
form factors |14| at both vertices 



t-lmd+v/ 2 2 2\ ^ g| g| {qf + g| + g 2 ) + P^{gf, gf, g\, _ , 



fe?, <ll Ql) = h {ql + q 2 2 ) 2 + /l 2 g 2 q\ + /i 3 (g? + q\) q\ + K q\ 

+h {ql + ql) + h ql + h 7 , (27) 

with g| = {qi + q^) 2 - In the spirit of resonance chiral theory [28J a Lagrangian with two 
multiplets of vector resonances was proposed recently in Ref. [JT] and references therein, 
which reproduces the above LMD+V Ansatz and which fulfills all the QCD short-distance 
constraints for the (VVP) Green's function. 

The constants hi in the Ansatz for ^md+v j n gq_ f|2(jp are determined as follows. 
The normalization with the pion decay amplitude tt° — > 77 in Eq. ffTUj) yields h 7 = 
-N c M^M^ 2 /{Att 2 F^) - h 6 ml - h 4 m 4 n = -14.83 GeV 6 - h % m\ - h 4 m*, where we used 
M Vl =M P = 775.49 MeV and My 2 = M p , = 1.465 GeV [32]. Note that in Refs. [HE] the 
small corrections proportional to the pion mass were dropped, assuming that the \hi\ are 
of order 1 — 10 in appropriate units of GeV. The Brodsky-Lepage behavior from Eq. ( fTTl) 
can be reproduced by choosing hi = GeV 2 . Furthermore, in Ref. [H] a fit to the CLEO 
data for the on-shell form factor ^ L M ? +V (m 2 , -Q 2 , 0) was performed, with the result 
/15 = (6.93 ± 0.26) GeV 4 — hj,m\. Again, the correction proportional to the pion mass 
was omitted in Refs. [HIE]. As pointed out in Ref. [21], the constant h 2 can be obtained 
from the higher-twist corrections in the OPE. Comparing with Eq. ( fTTl) yields the result 
h 2 = -4 {M\ + MyJ + (16/9) 5 2 ~ -10.63 GeV 2 . 

Within the LMD+V framework, the vector-tensor two-point function reads [H] 

cvx = (29) 

where we fixed the constant cvt using Eq. fl2T|) . As shown in Ref. [13] the OPE from 
Eq. (TTHD for leads to the relation 

/ii + /i 3 + /i 4 = 2c VT - (30) 
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As noted above, the value of the magnetic susceptibility xil 1 ) an d the relevant scale p are 
not precisely known. However, the LMD estimate ^ LMD = —2/My = —3.3 GeV -2 is close 
to xO = 1 GeV) = -(3. 15 ±0.30) GeV" 2 obtained in Ref. [50] using QCD sum rules with 
several vector resonances p,p', and p" . Assuming that the LMD/LMD+V framework is 
self-consistent, we will therefore take x — (—3.3 ±1.1) GeV~ 2 in our numerical evaluation, 
with a typical large- N c uncertainty of about 30%. This translates into the constraint 
h 3 + hi = (-4.3 ± 1.4) GeV 2 , corresponding to c VT = (-2.13 ± 0.71) GeV 2 . We will vary 
h 3 in the range ±10 GeV 2 and determine from Eq. flHUl) and vice versa. 

Note that using the off-shell LMD+V form factor at the external vertex leads to a 
short-distance behavior in the full light-by- light scattering contribution which at least 
qualitatively agrees with the OPE constraint derived in Ref. [21]. As stressed earlier, 
Ref. |21j only considers the pion-pole contribution with on-shell form factors; therefore 
a direct quantitative comparison with our approach is not possible. Nevertheless, taking 
first q\ ~ q\ ^> gf and then g| large, one obtains, together with the pion propagator in 
Eq. ([2]) [in the term with T2], an overall 1/gf behavior for the pion-exchange contribution, 
since at the external vertex we have |55] 



-r-LMD+V/,,2 „2 n\ q 3^°° h l + h 3 + h 4 _ 2c VT _ , , 

F/"°V7 W3>?3,UJ - M\M\ ~ M*Ml~ X - [6 } 

In the derivation we used Eqs. (I2"§]) and (I3"UI) . see also Eq. (122"]) . This l/q% behavior is 
as expected from Eq. (!25l) . reproduced earlier from Ref. [21]. On the other hand, if we 
would use a constant form factor proportional to the WZW term at the external vertex 
as proposed in Ref. [21] . we would get [55] 

where for simplicity we considered the chiral limit. That means that with the value of 
X = — Vc/(47r 2 F 2 ) from Ref. [13] we would in Eq. ( J3T1) precisely satisfy the short-distance 
constraint from Ref. [21] . 

The coefficient he in the LMD±V Ansatz is undetermined as well. We can obtain some 
indirect information on its size and sign in the following way. Low-energy constants in 
chiral Lagrangians can be estimated by starting with some resonance Lagrangian and then 
integrating out the heavy resonance states, usually at tree level. In particular for low- 
energy constants which are given by the exchanges of vector and axial-vector mesons, this 
procedure works in general quite well [561128]. Although for instance for vector mesons one 
can write down many different Lagrangians, it was shown in Ref. [28J that imposing QCD 
short- distance constraints on the resonance Lagrangian itself leads to unique estimates 
for the low-energy constants at order p 4 in the chiral Lagrangian. At order p 6 this is not 
true anymore [TB| IT4"]; nevertheless, it seems reasonable to reduce the model dependence 
by imposing again short-distance constraints on such resonance Lagrangians. 

Usually, only the exchange of the lightest resonance state in each channel is considered 
in this approach. One expects, however, some corrections to these estimates, with a typical 
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large- Nc error of about 30%, if also the exchanges of heavier resonance states are taken 
into account. In Ref. [H] this was shown to be true for one of two linear combinations 
of low-energy constants from the chiral Lagrangian of odd intrinsic parity at order p 6 
which enter in the low-energy expansion of the Green's function (VVP). With the LMD 
and LMD+V Ansatze for this Green's function, the relevant combination of low-energy 
constants is given by 

/i LMD _ F n N C _ 1 u 10 4 ,„o\ 

Ay * 2 ~ 8M 4 32n 2 M 2 - i - U Fl ' 



A lmd+v _ fV ri 5 N c f M 2 \ _ 10" 

V 'P 2 32tt 2 M 2 i \ M\ ) F* ' 



The constant h§ which enters A^P +V is directly related to the Brodsky-Lepage behavior 
of the form factor. Even though this falloff behavior cannot be reproduced with the LMD 
Ansatz, the change in the low-energy constant when going from LMD to LMD+V is only 
about 20%, well within the expected large- Nc uncertainty. 

On the other hand, the coefficient he determines a second linear combination of low- 
energy constants at order p 6 : 

F 2 10~ 4 F 2 

A^ q? = = -0.26 — , = ___zl_/, 6 . (35) 

Note that using the resonance Lagrangian of Ref. [H7] , one would obtain Ay S , +q -, 2 = in- 
stead. However, this resonance Lagrangian in general fails to reproduce the short-distance 
constraints from QCD, in contrast to the LMD and LMD+V Ansatze; see Ref. [H]. In 
particular, the prediction for Ay^ p+q y in the LMD model follows directly from the im- 
plementation of these short-distance constraints. Note that there is no problem with the 
short- distance behavior for the LMD form factor in the relevant channel where at low en- 
ergies A Vj (p+ q )2 enters. If we would assume a 30% error on the LMD estimate in Eq. ASS]) , 
we would obtain the quite narrow range he = M4(l ±0.3) = (4.6 ± 1.4) GeV 4 . However, 
this procedure might underestimate the potential variation of he, since the low-energy 
constant happens to be small compared to v4y^°; see Eq. ( 1331) . The magnitude 

of the shift of A v , p 2 when going from LMD to LMD+V is -0.25 (10^ 4 /F 2 ). That is, 
the shift is of the same size as ^v^ +q \2 itself. Assuming again that the LMD/LMD+V 
framework is self-consistent, but, to be conservative, allowing for a 100% uncertainty of 

A vS?+ q )*i we § et the ran S e ^6 = (5 ± 5) GeV 4 . 

Of course, the uncertainties of the values of the undetermined parameters h 3 , h$ and 
he and of the magnetic susceptibility x(aO is a drawback when using the off-shell LMD+V 
form factor and will limit the precision of the final estimate. 

The integral to be performed in Eq. (j2J) is eight-dimensional, thereof 3 integrations can 
be done trivially. In general, one then has to deal with a five-dimensional integration over 
3 angles and 2 moduli. We have performed these integrations numerically after a rotation 
to Euclidean momenta using the program VEGAS [58J . As a check we have reproduced the 
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values of a*; byL;7r for various form factors which have been used earlier in the literature. 

For instance, using a simple VMD form factor, we obtain o^vmd = 57 x 10" 11 for 
m M = 105.658369 MeV and m,o = 134.9766 MeV and with the value M p = 775.49 MeV. 

The results for a^ byL;7r for some selected values of h 3 , hi and he, varied in the ranges 
discussed above, for x — —3.3 GeV" 2 , hi = GeV 2 , h 2 = —10.63 GeV 2 and h 5 = 
6.93 GeV 4 - h 3 m 2 are collected in Table [TJ 

Table 1: Results for a^ jL ' n ° x 10 11 obtained with the off-shell LMD+V form factor for 
X = —3.3 GeV" 2 and the given values for h 3 , hi and he- When varying h 3 (upper half of 
the table), the parameter hi is fixed by the constraint in Eq. ( 1301) . In the lower half the 
procedure is reversed. The values of the other parameters are given in the text. 





h e = GeV 4 


h 6 = 5 GeV 4 


h 6 = 10 GeV 4 




= -10 GeV 2 


68.4 


74.1 


80.2 


h 


= GeV 2 


66.4 


71.9 


77.8 


h 


= 10 GeV 2 


64.4 


69.7 


75.4 


hi 


= -10 GeV 2 


65.3 


70.7 


76.4 


h^ 


= GeV 2 


67.3 


72.8 


78.8 


hi 


= 10 GeV 2 


69.2 


75.0 


81.2 



Varying x by ±1.1 GeV" 2 changes the result for a^ byL;7r ° by ±2.1 x 10" 11 at most. One 
observes from the table that the uncertainty in he affects the result by up to ±6.4 x 10~ u . 
If we would use instead h e = (0 ± 10) GeV 4 , the result would vary by about ±12 x 10" 11 
around the central value. The variation of aJ; byL;7r with h 3 [with /i 4 determined from 
the constraint in Eq. ( J30l) or vice versa] is much smaller, at most ±2.5 x 10" 11 . The 
variation of h$ by ±0.26 GeV 4 only leads to changes of ±0.6 x 10" 11 in the final result. 
Within the scanned region, we obtain a minimal value of a^ byL;?r = 63.2 x 10 -11 for x — 
-2.2 GeV" 2 , h 3 = 10 GeV 2 , and h 6 = GeV 4 and a maximum of aJ; byL ^ = 83.3 x 10" 11 for 
X = —AA GeV -2 , hi = 10 GeV 2 , and he = 10 GeV 4 . In the absence of more information 
on the precise values of the constants h 3 ,hi and h 6 , we take the average of the results 
obtained with h 6 = 5 GeV 4 for h 3 = GeV 2 , i.e. 71.9 x 10" 11 , and for h 4 = GeV 2 , 
i.e. 72.8 x 10" 11 , as our central value, 72.3 x 10" 11 . To estimate the error, we add all 
the uncertainties from the variations of x> h 3 (or hi), h$ and he linearly to cover the full 
range of values obtained with our scan of parameters. Note that the uncertainties of x 
and the coefficients h 3 , hi and he do not follow a Gaussian distribution. In this way we 
obtain our final estimate 

a LbyL;^ = ^ ± ^ x ^-11 ^ 

We think the 16% error should fairly well describe the inherent model uncertainty using 
the off-shell LMD+V form factor. In order to facilitate future updates of our result in 
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case some of the parameters hi in the LMD+V Ansatz in Eq. (1261) or the value (and 
the relevant scale) of the magnetic susceptibility xil 1 ) wm be known more precisely, we 
present in the appendix a parametrization of a^ byL;7r for arbitrary coefficients hi. 

As far as the contribution to a M from the exchanges of the other light pseudoscalars 
77 and 7]' is concerned, it is not so straightforward to apply the above analysis within 
the LMD+V framework to these resonances. In particular, the short- distance analysis in 
Ref. [14J was performed in the chiral limit and assumed octet symmetry. For the rj the 
effect of nonzero quark masses has definitely to be taken into account. Furthermore, the 
rj' has a large admixture from the singlet state and the gluonic contribution to the axial 
anomaly will play an important role. We therefore resort to a simplified approach which 
was also adopted in other works [6l [7J [8j [H [21] and take a simple VMD form factor 

.77VMD / 2 2 2\ _ N C My- My _ , 

U7l 2 Fps 52 _ M 2 ) {q 2 _ M 2 j > ^~V,V, W 

normalized to the experimental decay width T(PS — > 77). We can fix the normalization 
by adjusting the (effective) pseudoscalar decay constant Fps in Eq. ( 1371) . Using the latest 
values T(r] -> 77) = (0.510 ± 0.026) keV and r(r/ -> 77) = (4.30 ± 0.15) keV from 
Ref. [32], one obtains F VfiS = 93.0 MeV with m v = 547.853 MeV and = 74.0 MeV 

with mrf = 957.66 MeV. However, we do not follow the approach of Ref. [21] and will 
also take a VMD form factor at the external vertex. 

Note that the on- and off-shell VMD form factors are identical, since the form factor 
does not depend on the momentum g| which flows through the pion leg. The prob- 
lem with the VMD form factor is that the damping is too strong as it behaves like 
Jv) 7 » 7 » (m 2 , —Q 2 , —Q 2 ) ~ 1/Q 4 , instead of ~ l/Q 2 deduced from the OPE, see Eq. ffTBj) . 
This effect might lead to an underestimating of the contribution. However, the relevant 
integrals for a^ byL;PS do not seem to be very sensitive to the correct asymptotic behavior 
for large momenta. This can be seen from the weight functions which multiply the form 
factors in the integral and which are displayed in Ref. [9]. It seems more important to have 
a good description at small and intermediate energies below 1 GeV, e.g. by reproducing 
the slope of the form factor JF PS7 » 7 (— Q 2 , 0) at the origin. The CLEO Collaboration [3"4] 
has made a fit of the on-shell form factors T^^i—Q 2 , 0) and ^' 7 * 7 (— Q 2 , 0), normalized 
to the corresponding experimental width T(PS — > 77), using a VMD Ansatz with an 
adjustable parameter Aps in place of the vector- meson mass My in Eq. (j3~Tj) . Taking 
their values = (774 ± 29) MeV and A„/ = (859 ± 28) MeV, we then obtain the results 
a Lb y L;r; _ ^ 5 x ^g-ii and a Lb y L ; r/ _ ^.5 x 10" 11 , which update the values given in 

Ref. [9jO Only a more detailed analysis, along the line of the LMD+V framework, will 
show whether these values are realistic. Thus, adding up the contributions from all the 
light pseudoscalar exchanges (ir°, 77, 77'), we obtain the estimate 

a LbyL;PS = (gg ± ^ x ^-1^ (3g) 

6 If we use a constant (WZW) form factor at the external vertex, as proposed in Ref. [3T], we would 
obtain a ^y h ^-P° lc = 21.5 x 1CT 11 and a LbyL ;l /- P oic = 20 .i x 10" 11 . Note that these values are somewhat 
larger than a^y^-P 010 = 18 x 10" 11 and a LbyL;V- P oic = 18 x 10" 11 given in Ref. |2T] . 
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where we have assumed a 16% error, as inferred above for pion-exchange contribution 
using the off-shell LMD+V form factorQ 

5 Discussion and conclusions 

Following the observation in Refs. [251 EE] we have reevaluated the pion-exchange contribu- 
tion to hadronic light-by-light scattering in the muon g — 2 using fully off-shell form factors 
JF^o*^. at both vertices. We used a model based on the large- Nc QCD framework with 
two multiplets of vector mesons (LMD+V) [H] which fulfills all QCD short-distance con- 
straints on the form factor and reproduces the experimentally confirmed Brodsky-Lepage 
behavior. We also derived a new short- distance constraint on the form factor at the ex- 
ternal vertex, relating it to the quark condensate magnetic susceptibility \. The obtained 
value a^ byL;7r ° = (72 ± 12) x 1CT 11 replaces the result obtained in Ref. [9] with on-shell 
LMD+V form factors at both vertices. Adding the contribution from the exchanges of the 
i] and rf evaluated with simple VMD form factors, we obtain a^ byL;PS = (99 ± 16) x 10~ n 
for the sum of all light pseudoscalars. These values for the pion and all pseudoscalars are 
about 20% larger than the estimates obtained in Refs. [13 LH E] which used other hadronic 
models for the form factor. 

As mentioned earlier, the identification of individual contributions, like pion exchange, 
in hadronic light-by- light scattering is model-dependent as soon as one uses off-shell form 
factors. We view our evaluation as being a part of a full calculation based on a resonance 
Lagrangian, which fulfills all the relevant QCD short-distance constraints, along the lines 
of the resonance chiral theory approach developed in Ref. |28j . 

We would like to stress that although our result for the pion-exchange contribution is 
not too far from the value a Lb y L;7r - P oie _ (75.5 + 6.7) x 10~ n given in Ref. [21] JU this is pure 
coincidence. We have used off-shell LMD+V form factors at both vertices, whereas the 
authors of Ref. [21] evaluated the pion-pole contribution using the on-shell LMD+V form 
factor at the internal vertex and a constant (WZW) form factor at the external vertex. On 
the other hand, as has been observed in Refs. [6j U\ El El [22] , it is the region of momenta 
below about 2 GeV which gives the bulk of the contribution to the final result in the pion- 
exchange or pion-pole correction to hadronic light-by- light scattering. This is also clearly 
visible from the weight functions that multiply the form factors in the integrals and which 
have been presented in Ref. [§]; see also Ref. [22]- They have a peak around 0.5 GeV. 

Applying the same procedure to the electron, we obtain = (2.98 ± 0.34) x 1(T 14 with off- 

shell LMD+V form factors at both vertices. This number supersedes the value given in Ref. [9]. Note 
that the naive rescaling aj byL;T (rescaled) = (■m e /m fl ) 2 a ! ^ byL '' 7r = 1.7 X 10~ 14 yields a value which is 
almost a factor of 2 too small. Our estimates for the other pseudoscalars contributions using VMD 
form factors at both vertices are a^ byh ^ = 0.49 x 10~ 14 and <+ byL;,? ' = 0.39 x 10" 14 . Therefore we 
get a^ byL;PS = (3.9 ± 0.5) x 10~ 14 , where the relative error of about 12% is again taken over from the 
pion-exchange contribution. 

8 Actually, using the on-shell LMD+V form factor at the internal vertex with /12 = —10 GeV 2 and 
/15 = 6.93 GeV 4 and a constant (WZW) form factor at the external vertex, we obtain 79.8 x 10~ n , close 
to the value 79.6 x lO" 11 given in Ref. [22] and 79.7 x 10~ n in Ref. [29]. 
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Therefore, as long as the absolute values of the model parameters and h e , which 

control the off-shellness of the pion in the LMD+V form factor in Eq. fl2U|) J^I are not too 
large, i.e. below about 10 in appropriate units of GeV, one obtains a result which will not 
be too far from the one obtained with on-shell LMD+V form factors. We have given some 
arguments for our choice of the parameters hi and the ranges in which we vary them and 
they fulfill this constraint on their size. Recall that the constant term in the numerator 
of the form factor in Eq. fT2U|) has the value h 7 ~ —14.8 GeV 6 . As pointed out before, our 
Ansatz for the neutral pion contribution to hadronic light-by-light scattering with two off- 
shell LMD+V form factors agrees qualitatively with the short-distance behavior derived 
in Ref. [2T] . However, since only the pion-pole contribution was considered throughout 
that paper, a direct quantitative comparison is not possible. 

Recently, an evaluation of the pion-exchange contribution appeared [2H] which uses a 
nonlocal chiral quark model for the off-shell form factor jF 7r0 * 7 * 7 *. In that model, off-shell 
effects of the pion always lead to a rather strong damping in the form factor and the result 
a Lb y L;7T _ + 2) x 10~ n is therefore smaller than the pion-pole contribution obtained 
in Ref. [21]. Although we also get a value which is (slightly) smaller than the pion-pole 
contribution, our result depends on the chosen model parameters, i.e. the constants hi 
in the LMD+V Ansatz and on the value for the magnetic susceptibility xif 1 )- We have 
given arguments for our preferred choice of these parameters. In some other corner of 
the parameter space one can, however, obtain a result which is larger than the pion-pole 
contribution, i.e. we get a maximal value of 

Q LbyL;^ = gg 3 x 1Q -11 ^ thg scanne d 

range. Of 

course, any additional information to pin down these model parameters and the "correct" 
value of xif 1 ) an d the relevant scale \i would be highly welcome to obtain a more precise 
prediction for the pion-exchange contribution. At this point it is not clear whether the 
nonlocal chiral quark model used in Ref. [29] or our LMD+V model for the form factor 
better represents the strongly interacting region of QCD below about 2 GeV. At least 
the LMD+V form factor fulfills all the relevant QCD short-distance constraints. In any 
case, we think that the error of ±2 x 10~ n given in Ref. [2H] probably underestimates the 
inherent uncertainly of any hadronic model. 

In Ref. [21] an improved evaluation of the axial-vector contribution to hadronic light- 
by-light scattering was given compared to Refs. [HI [3, [8], with the result (22 ± 5) x 1CT 11 . 
Note, however, that this seems to be again only the pole contribution. Furthermore, 
Ref. [0] obtained the following results for the remaining contributions: (—7 + 2) x 10~ n for 
scalar exchanges, (—19+13) x 10~ n for the dressed pion and kaon loops and (21+3) x 10 -11 
for the dressed quark loops. These estimates have more conservative errors than those 
given in Refs. pD, [8]. Furthermore, the scalar-exchange contribution is not evaluated in 
the latter references. If we combine our value for the pseudoscalars with these results, we 
obtain the new estimate 

a LbyL;had = ( n g ± ^ x jq-U (gg) 

9 Note, however, that even if /13, /14 and are put to zero, which actually violates the short-distance 
constraint from Eq. (130|) . one does not recover the on-shell LMD+V form factor because of a term 
proportional to q\q\{q\ + 92) 2 in the numerator in Eq. (|26| . 
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for the total hadronic light-by-light scattering contribution to the anomalous magnetic 
moment of the muon. To be conservative, we have added all the errors linearly, as has 
become customary in recent years, and rounded up the obtained value ±39 x 10 -11 . In 
the very recent review [21] the central values of some of the individual contributions to 
hadronic light-by- light scattering are adjusted and some errors are enlarged to cover the 
results obtained by various groups which used different models. The errors are finally 
added in quadrature to yield the estimate a L b y L ; had — (105 ± 26) x 10~ u . Note that the 
dressed light quark loops are not included as a separate contribution in Ref. [21]. They 
are assumed to be already covered by using the short-distance constraint from Ref. [21] 
on the pseudoscalar-pole contribution. 

Some progress has been achieved in recent years to better understand the hadronic 
light-by-light scattering contribution to the muon g — 2. We hope that our new short- 
distance constraint on the off-shell form factor jF 7r0 * 7 » 7 (g 2 , q 2 , 0) at the external vertex will 
further help to control the numerically dominant pion-exchange contribution. We should 
not forget, however, that the contribution of the exchanges of rj and rf are theoretically 
not that well understood. We have simply used VMD form factors as has been done in 
most other evaluations. A new analysis, along the lines of the approach for the pion, is 
definitely needed. Finally, as stressed in Refs. [2"2"| ]24j a better control of the numerically 
subdominant but non-negligible other contributions is also needed, if we fully want to 
profit from a potential future g — 2 experiment. 
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Appendix 

We provide in this appendix a parametrization of a^ byL;7r ° for arbitrary coefficients hi in 
the LMD+V Ansatz for the off-shell form factor jF 7r0 * 7 » 7 . in Eq. (|2"B"]) . This will facilitate 
future updates of our result for the pion-exchange contribution in Eq. ( l36|) in case some 
of the parameters hi or the value (and the relevant scale) of the magnetic susceptibility 
x(/x) will be known more precisely. 

Measuring the parameters hi in appropriate units of GeV, i.e. defining hi = hi/GeV 2 
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for z 



1, 2, 3, 4, hi = hi/GeV 4 for i = 5, 6 and /17 = h 7 /GeV 6 , we can write 



LbyL;7r' 



(Al) 



_i=l j=l j=j 



where the coefficients Cj and c^- are given in Table [2] 



Table 2: Values of the coefficients q and Cy which appear in Eq. (1A11) . 





% = 1 


% = 2 


z = 3 


z = 4 


z = 5 


z = 6 


z = 7 


Cj X 10 4 


-1.4530 





-1.4530 


-1.4530 


0.4547 


0.4547 


-1.2048 


Cy X 10 4 


J = l 


j = 2 


J = 3 


i = 4 


J = 5 


j = 6 


J = 7 


I = 1 


0.4447 


0.0729 


0.7428 


0.7120 


-0.3620 


-0.3332 


0.9916 


z = 2 







0.0730 


0.0729 


-0.0557 


-0.0557 


0.2221 


% = 3 






0.2980 


0.5653 


-0.1967 


-0.1679 


0.1796 


i = 4 








0.2672 


-0.1679 


-0.1391 


0.1162 


z = 5 










0.1215 


0.1796 


-0.8072 


z = 6 












0.0581 


-0.3052 


i = 7 














1.6122 



This representation follows immediately from the general expression for the pion- 
exchange contribution in Eq. (j2J) and from the LMD+V Ansatz for the form fac- 
tor in Eq. ( !26|) . Because of our new short-distance constraint at the external ver- 
tex, the parameters /ii, /13 and are not independent, but must obey the relation 
hi + h 3 + hi = My x My 2 x', see Eq. ( 1301) . Note the absence of a term without the con- 
stants hi in Eq. (IA10 . This follows from the fact that at the external vertex with the 
soft photon the form factor ^ r 7r o* 7 * 7 (q , |, qf, 0) enters, e.g. in the term with T 2 in Eq. ([2]). 
This also leads to c 2 = and c 22 = in Table [21 For evaluating the integrals, we used 
ttv = 105.658369 MeV, m^o = 134.9766 MeV, F n = 92.4 MeV, M Vl = M p = 775.49 MeV 
and My 2 = M p i = 1.465 GeV. The decay constant F w only enters as an overall factor in 



Note, however, that some of the model parameters hi are quite well determined from 
experimental data and theoretical constraints. The normalization with the pion decay 



amplitude vr° -> 77 yields h 7 = -N c M Vi M v J(A-n 2 F%) - h^m\ - h±m%. The Brodsky- 
Lepage behavior can be reproduced by choosing h\ = GeV 2 . Furthermore, a fit to 
the CLEO data for the on-shell form factor j^° +v (m;, -Q 2 ,0) leads to h b = (6.93 ± 
0.26) GeV 4 — h 3 m 2 . Finally, the constant h 2 can be obtained from higher-twist corrections 
in the OPE, with the result h 2 = -4 (M 2 x + M 2 J + (16/9) 5 2 ~ -10.63 GeV 2 . 
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If we use the above informations to fix hi,h 2 ,h 5 and h 7 , we obtain the simplified 
expression 



503.3764 - 6.5223 h 3 - 5.0962 h A + 7.8557 h 6 
+0.3017^ + 0.5683 h 3 h 4 - 0.1747^ h 



-0.2672 h\ - 0.1411 h A h e + 0.0642 h\ 



x 10" 4 , (A2) 



where only h 3 , /i 4 and Hq enter as free parameters. Note again, however, that h 3 and h A 
are not independent, but now obey the relation h 3 + h A = My M V2 x- 
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